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We discuss the results of laboratory measurements and theoretical models concerning the 
aggregation of dust in protoplanetary disks, as the initial step toward planet formation. Small 
particles easily stick when they collide and form aggregates with an open, often fractal structure, 
depending on the growth process. Larger particles are still expected to grow at collision 
velocities of about Im/s. Experiments also show that, after an intermezzo of destructive 
velocities, high collision velocities above lOm/s on porous materials again lead to net growth of 
the target. Considerations of dust-gas interactions show that collision velocities for particles not 
too different in surface-to-mass ratio remain limited up to sizes about Im, and growth seems 
to be guaranteed to reach these sizes quickly and easily. For meter sizes, coupling to nebula 
turbulence makes destructive processes more likely. Global aggregation models show that in 
a turbulent nebula, small particles are swept up too fast to be consistent with observations of 
disks. An extended phase may therefore exist in the nebula during which the small particle 
component is kept alive through collisions driven by turbulence which frustrates growth to 
planetesimals until conditions are more favorable for one or more reasons. 



1. INTRODUCTION 

The growth of dust particles by aggregation stands at 
the beginning of planet formation. Whether planetesimals 
form by incremental aggregation, or through gravitational 
instabilities in a dusty sublayer, particles have to grow and 
settle to the midplane regardless. On the most basic level, 
the physics of such growth is simple: Particles collide be- 
cause relative velocities are induced by random and (size- 
dependent) systematic motions of grains and aggregates in 
the gaseous nebula surrounding a forming star. The details 
are, however, highly complex. The physical state of the 
disk, in particular the presence or absence of turbulent mo- 
tions, set the boundary conditions. When particles collide 
with low velocities, they stick by mutual attractive forces, 
be it simple van der Waals attraction or stronger forces 
(molecular dipole interaction in polar ices, or grain-scale 
long-range forces due to charges or magnetic fields). While 
the lowest velocities create particle shapes governed by the 
motions alone, larger velocities contribute to shaping the 
aggregates by restructuring and destruction. The ability to 
internally dissipate energy is critical in the growth through 
intermediate pebble and boulder sizes. In this review we 
will concentrate on the physical properties and growth char- 
acteristics of these small and intermediate sizes, but also 



make some comments on the formation of planetesimals. 

Relative velocities between grains in a protoplanetary 
disk can be caused by a variety of processes. For the small- 
est grains, these are dominated by Brownian motions, that 
provide relative velocities in the mm/s to cm/s range for 
(sub)micron sized grains. Larger grains show systematic 
velocities in the nebula because they decouple from the gas, 
settle vertically, and drift radially. At 1 AU in a solar nebula, 
these settling velocities reach m/s for cm-sized grains. Ra- 
dial drift becomes important for even larger particles and 
reaches lO's of m/s for m-sized bodies. Finally, turbu- 
lent gas motions can induce relative motions between par- 
ticles. For details see for example Weidenschilling (1977; 
1984), Weidenschilling and Cuzzi (1993), Cuzzi and Hogan 
(2003). 

The timescales of growth processes and the density and 
strength of aggregates formed by them, will depend on the 
structure of the aggregates. A factor of oveiTiding im- 
portance for dust-gas interactions (and therefore for the 
timescales and physics of aggregation), for the stability of 
aggregates, and for optical properties alike is the structure 
of aggregates as they form through the different processes. 

The interaction of particles with the nebula gas is deter- 
mined primarily by their gas drag stopping time tg which is 
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Fig. 1. — Projected area of aggregates as a function of aggregate 
size and fractal dimension, normalized to the cross section of a 
compact particle with the same mass. 



given by 
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where m is the mass of a particle, v its velocity relative to 
the gas, a the average projected surface area, pg is the gas 
density, Cg is the sound speed, and Fi^c is the drag force. 
The second equal sign in eq.(l) holds under the assumption 
that particles move at sub- sonic velocities and that the mean 
free path of a gas molecule is large compared to the size of 
the particle (Epstein regime). In this case, the stopping time 
is proportional to the ratio of mass and cross section of the 
particle. For spherical non-fractal (i.e. compact or porous) 
particles of radius o and mass density ps. this can be written 
as = aps/cpg. Fractal particles are characterized by the 
fact that the average density of a particle depends on size 
in a powerlaw fashion, with a power (the fractal dimension 
Di) smaller than 3. 
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For large aggregates, this value can in principle be mea- 
sured for individual particles. For small particles, it is often 
more convenient to measure it using sizes and masses of a 
distribution of particles. 

Fractal particles generally have large surface-to-mass ra- 
tios; in the limiting case of long linear chains (Df = 1) of 
grains with radii ao, cr/m approaches the constant value 
Stt/ (16aoPs). This value differs from the value for a single 
grain 3/(4aops) by just a factor 7r/4. Fig. 1 shows how the 
cross section of particles varies with their mass for differ- 
ent fractal dimensions. It shows that for aggregates made of 
10000 monomers, the surface-to-mass ratio can easily differ 
by a factor of 10. An aggregate made from O.l/um particles 
with a mass equivalent to a lO/Lim particle consists of 10^ 
monomers and the stopping time could vary by a factor of 
order 100. Just how far the fractal growth of aggregates 
proceeds is really not yet known. 

This review is organized as follows: In section 2 we 
cover the experiments and theory describing the basic 



growth processes of dust aggregates. In section 3 we 
discuss particle-gas interactions and the implications for 
inter-particle collision velocities as well as planetesimal 
formation. In section 4 we describe recent advances in the 
modeling of dust aggregation in protoplanetary disks and 
observable consequences. 

2. DUST AGGREGATION EXPERIMENTS AND 
THEORY 

2.1. Interactions between individual dust grains 

2.1.1 Interparticle adhesion forces. Let us assume that the 
dust grains are spherical in shape and that they are elec- 
trically neutral and non-magnetic. In that case, two grains 
with radii ai and 02 will always experience a short-range at- 
traction due to induced dielectric forces, e.g. van der Waals 
interaction. This attractive force results in an elastic defor- 
mation leading to a flattening of the grains in the contact 
region. An equilibrium is reached when the attractive force 
equals the elastic repulsion force. For smidl, hard grains 
with low surface forces, the equilibrium contact force is 
given by (Derjaguin et ai, 1975) 
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where 7s and R denote the specific surface energy of the 
grain material and the local radius of surface curvature, 
given hy R = a\a2/{a\ + 02), respectively. Measurements 
of the separation force between pairs of Si02 spheres with 
radii a between 0.5 jim and 2.5 /xm (corresponding to re- 
duced radii R = 0.35 .. . 1.3 /um) confirm the vahdity of 
Eq. 3 (Heimetal., 1999). 

2.7.2 Interparticle rolling-friction forces. Possibly the most 
important parameter influencing the structure of aggregates 
resulting from low velocity colUsions is the resistance to 
rolhng motion. If this resistance is very strong, both aggre- 
gate compaction and internal energy dissipation in aggre- 
gates would be very difficult. Resistance to rolling first of 
aU depends strongly on the geometry of the grains. If grains 
contain extended flat surfaces, contact made on such loca- 
tions could not be moved by rolling - any attempt to roll 
them would inevitably lead to breaking the contact. In the 
contact between round surfaces, resistance to roUing must 
come from an asymmetric distribution of the stresses in the 
contact area. Without external forces, the net torque exerted 
on the grains should be zero. Dominik and Helens (1995) 
showed that the pressure distribution becomes asymmetric, 
when the contact area is slightly shifted with respect to the 
axis cormecting the curvature centers of the surfaces in con- 
tact. The resulting torque is 
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where Ocontact.o is the equihbrium contact radius, Ocontact 
the actual contact radius due to pressure in the vertical di- 
rection, and ^ is the displacement of the contact area due to 
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the torque. In this picture, energy dissipation, and therefore 
friction, occurs when the contact area suddenly readjusts af- 
ter it has been displaced because of external forces acting on 
the grains. The friction force is proportional to the pull-off 

force Fc- 

Heim et al. (1999) observed the reaction of a chain 
of dust grains using a long-distance microscope and mea- 
sured the applied force with an Atomic Force Microscope 
(AFM). The derived rolling-friction forces between two 
Si02 spheres with radii of o = 0.95 are Froii = 
(8.5 ± 1.6) • 10~^° N. If we recall that there are two grains 
involved in rolling, we get for the rolling-friction energy, 
defined through a displacement of an angle 7r/2 

£;roii = ^aFrou = O(10-i5j) . (5) 

Recently, the rolling of particle chains has been observed 
under the scanning electron microscope while the contact 
forces were measured simultaneously {Heim et al., 2005). 

2.1.3 Sticking efficiency in single grain collisions. The dy- 
namical interaction between small dust grains was derived 
by Poppe et al. (2000a) in an experiment in which single, 
micrometer-sized dust grains impacted smooth targets at 
various velocities (0 . . . 100 m/s) under vacuum conditions. 
For spherical grains, a sharp transition from sticking with 
an efficiency of /? w 1 to bouncing (i.e. a sticking effi- 
ciency of /? = 0) was observed. This threshold velocity 
is Vs ~ 1.2 m/s for a = 0.6/zm and Vg ~ 1.9 m/s for 
a = 0.25 /im. It decreases with increasing grain size. The 
target materials were either polished quartz or atomically- 
smooth (surface-oxidized) siUcon. Currently, no theoreti- 
cal explanation is available for the threshold velocity for 
sticking. Earher attempts to model the low-velocity impact 
behavior of spherical grains predicted much lower sticking 
velocities (Chokshi et al., 1993). These models are based 
upon impact experiments with "softer" polystyrene grains 
(Dahneke, 1975). The main difference becomes visible 
when studying the behavior of the rebound grains in non- 
sticking colhsions. In the experiments by Dahneke (1975) 
and also in those by Bridges et al. (1996) using macro- 
scopic ice grains, the behavior of grains after a bouncing 
colUsion was a unique function of the impact velocity, with 
a coefficient of restitution (rebound velocity divided by im- 
pact velocity) always close to unity and increasing mono- 
tonically above the threshold velocity for sticking. For 
harder, still spherical, Si02 grains {Poppe et al, 2000a), 
the average coefficient of restitution decreases considerably 
with increasing impact velocity. In addition to that, individ- 
ual grain impacts show considerable scatter in the coeffi- 
cient of restitution. 

The impact behavior of irregular dust grains is more 
complex. Irregular grains of various sizes and composi- 
tions show an overall decrease in the sticking probability 
with increasing impact velocity. The transition from /? = 1 
to /3 = 0, however, is very broad so that even impacts as 
fast as w 100 m/s can lead to sticking with a moderate 
probabihty. 



2.2. Dust aggregation and restructuring 

2.2.7 Laboratory and microgravity aggregation experi- 
ments. In recent years, a number of laboratory and mi- 
crogravity experiments have been carried out to derive the 
aggregation behavior of dust under conditions of young 
planetary systems. To be able to compare the experimen- 
tal results to theoretical predictions and to allow numer- 
ical modelling of growth phases that are not accessible 
to experimental investigation, "ideal" systems were stud- 
ied, in which the dust grains were monodisperse (i.e. all 
of the same size) and initially non-aggregated. Whenever 
the mean collision velocity between the dust grains or ag- 
gregates is much smaller than the sticking threshold (see 
section 2.1.3), the aggregates formed in the experiments 
are "fractal", i.e. £>f < 3 {Wurm and Blum, 1998; Blum 
et al, 1999; Blum et al, 2000; Krause and Blum, 2004). 
The precise value of the fractal dimension depends on the 
specific aggregation process and can reach values as low as 
Df = 1.4 for Brownian-motion driven aggregation {Blum 
et al, 2000; Krause and Blum, 2004; Paszun and Dominik, 
2006), Df = 1.9 for aggregation in a turbulent gas {Wurm 
and Blum, 1998), or Di = 1.8 for aggregation by gravi- 
tationally driven sedimentation in gas {Blum et al, 1999). 
It is inherent to a dust aggregation process in which aggre- 
gates with low fractal dimensions are formed that the mass 
distribution function is rather narrow (quasi-monodisperse) 
at any given time. In all reaUstic cases, the mean aggregate 
mass in follows either a power law with time t, i.e. fh oc f 
with 7 > {Krause and Blum, 2004) or grows exponen- 
tially fast, TO oc exp((5t) with (5 > {Wurm and Blum, 
1998) which can be verified in dust-aggregation models 
(see Section 2.2.2). 

As predicted by Dominik and Tielens (1995, 1996, 
1997), experiments have shown that at collision velocities 
near the velocity threshold for sticking (of the individual 
dust grains), a new phenomenon occurs {Blum and Wurm, 
2000). Whereas at low impact speeds, the aggregates' struc- 
tures are preserved in collisions (the so-called "hit-and- 
stick" behavior), the forming aggregates are compacted at 
higher velocities. In even more energetic collisions, the 
aggregates fragment so that no net growth is observable. 
The different stages of compaction and fragmentation are 
depicted in Fig. 2. 

2.2.2 Modelling of dust aggregation. The evolution of grain 
morphologies and masses for a system of initially monodis- 
perse spherical grains that are subjected to Brownian mo- 
tion has been studied numerically by Kempf et al (1999). 
The mean aggregate mass increases with time following a 
power law (see Section 2.2.1). The aggregates have frac- 
tal structures with a mean fractal dimension of Df = 1.8. 
Analogous experiments by Blum et al (2000) and Krause 
and Blum (2004), however, found that the mean fractal di- 
mension was £>£ = 1.4. Recent numerical work by ftjiZMn 
and Dominik (2006) showed that this lower value is caused 
by Brownian rotation (neglected by Kempf et al (1999)). 
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Fig. 2. — Dominating processes and associated energies in aggre- 
gate collisions, after Dominik and Tielens (1997) and Wurm and 
Blum (2000). i?br is the energy needed to break a contact, -Eroii 
is the energy to roll two grains through an angle 7r/2 and n is the 
total number of contacts in the aggregates. 



More chain-like dust aggregates can form if the mean free 
path of the colliding aggregates becomes smaller than their 
size, i.e. if the assumption of balUstic collisions breaks 
down and a random-walk must be considered for the ap- 
proach of the particles. The fractal dimension of thermally 
aggregating dust grains is therefore dependent on gas pres- 
sure and reaches an asymptotic value of Df = 1.5 for the 
low density conditions prevailing through most of the preso- 
lar nebula. Only in the innermost regions the densities are 
high enough to cause deviations. 

The experimental work reviewed in Section 2.2.1 can 
be used to test the applicability of theoretical dust aggre- 
gation models. Most commonly, the mean-field approach 
by Smoluchowski (1916) is used for the description of the 
number density n(m, t) of dust aggregates with mass m as 
a function of time t. Smoluchowski' s rate equation reads in 
the integral form 

dn{m,t) ^ r T^i ' '^ 

= K{m,m-m) (6) 

■n{m', t) n{m — m', t) dm' 

—n{m,t) / K {m' , m) n{m' , t) dm' . 
Jo 

Here, K{mi,m2) is the reaction kernel for aggregation of 
the coagulation equation 6. The first term on the rhs. of 
Eq. 6 describes the rate of sticking collisions between dust 
particles of masses m.' and m — m' whose combined masses 
are m (gain in number density for the mass m). The second 
term denotes a loss in the number density for the mass m 
due to sticking collisions between particles of mass m and 
mass m'. The factor 1/2 in the first term accounts for the 
fact that each pair colUsion is counted twice in the integral. 
In most astrophysical applications the gas densities are so 
low that dust aggregates collide ballistically. In that case. 



the kernel in Eq. 6 is given by 

K{mi,m2) = /3(mi, 7712; w) w(mi, 7712) CT(mi, TO2) , (7) 

where f3{mi,m2',v), u(mi,m2), and a{m,i,m2) are the 
sticking probability, the collision velocity, and the cross 
section for collisions between aggregates of masses m\ and 
m2, respectively. 

A comparison between numerical predictions from Eq. 
6 and experimental results on dust aggregation was given 
by Wurm and Blum (1998) who investigated dust aggrega- 
tion in rarefied, turbulent gas. Good agreement for both the 
mass distribution functions and the temporal behavior of the 
mean mass was found when using a sticking probability of 
P{mi,m2;v) = 1, a mass-independent relative velocity be- 
tween the dust aggregates and the expression by Ossenkopf 
(1993) for the collision cross section of fractal dust aggre- 
gates. Blum (2006) showed that the mass distribution of the 
fractal aggregates observed by Krause and Blum (2004) for 
Brownian-motion driven aggregation can also be modelled 
in the transition regime between free-molecular and hydro- 
dynamic gas flow. 

Analogous to the experimental findings for the colU- 
sional behavior of fractal dust aggregates with increasing 
impact energy (Blum and Wurm (2000), see Section 2.2.1), 
Dominik and Tielens (1997) showed in numerical experi- 
ments on aggregate collisions that with increasing coUision 
velocity the foUowing phases can be distinguished: hit-and- 
stick behavior, compaction, loss of monomer grains, and 
complete fragmentation (see Fig. 2). They also showed that 
the outcome of a collision depends on the impact energy, the 
rolling-friction energy (see Eq. 5 in Section 2.1.2) and the 
energy for the breakup of single interparticle contacts (see 
Section 2.1.1). The model by Dominik and Tielens (1997) 
was quantitatively confirmed by the experiments of Blum 
and Wurm (2000) (see Fig. 2). 

To analyze observations of protoplanetary disks and 
model the radiative transfer therein, the optical properties 
of particles are important (McCabe et ai, 2003; Ueta and 
Meixner, 2003; Wolf, 2003). Especially for particle sizes 
comparable to the wavelength of the radiation, the shape 
and morphology of a particle are of major influence for the 
way the particle interacts with the radiation. With respect 
to this, it is important to know how dust evolution changes 
the morphology of a particle. As seen above, in most cases 
dust particles are not individual monolithic solids but rather 
aggregates of primary dust grains. Numerous measure- 
ments and calculations have been carried out on aggregates 
(e.g. Kozasa et ai, 1992; Henning and Stognienko, 1996; 
Wurm and Schnaiter, 2002; Gustafson and Kolokolova, 
1999; Wurm et al, 2004a; Min et al, 2005). No simple 
view can be given within the frame of this paper. However, 
it is clear that the morphology and size of the aggregates 
will strongly influence the optical properties. 

2.2.3 Aggregation with long-range forces. Long range 
forces may play a role in the aggregation process, if grains 
are either electrically charged or magnetic. Small iron 



grains may become spontaneously magnetic if they are sin- 
gle domain {Nuth et ai, 1994; Nuth and Wilkinson, 1995), 
typically at sizes of a few tens of nanometers. Larger 
grains containing ferromagnetic components can be magne- 
tized by an impulse magnetic field generated during a light- 
ning discharge {Tiinyi et ai, 2003). For such magnetized 
grains, the colUsional cross section is strongly enhanced 
compared to the geometrical cross section {Dominik and 
NUbold, 2002). Aggregates formed from magnetic grains 
remain strong magnetic dipoles, if the growth process keeps 
the grain dipoles aligned in the aggregate (Niibold and 
Glassmeier, 2000). Laboratory experiments show the spon- 
taneous formation of elongated, almost Unear aggregates, 
in particular in the presence of an external magnetic field 
{Niibold et at., 2003). The relevance of magnetic grains to 
the formation of macroscopic dust aggregates is, however, 
unclear. 

Electric charges can be introduced through tribo-electric 
effects in collisions, through which electrons and/or ions 
are exchanged between the particles (Poppe et ai, 2000b; 
Poppe and Schrapler, 2005; Desch and Cuzzi, 2000). The 
number of separated elementary charges in a collision be- 
tween a dust particle and a solid target with impact energy 
i^c can be expressed by {Poppe et ai, 2000b; Poppe and 
Schrapler, 2005) 

The cumulative effect of many non-sticking collisions can 
lead to an accumulation of charges and to the build-up 
of strong electrical fields at the surface of a larger aggre- 
gate. In this way, impact charging could lead to electro- 
static trapping of the impinging dust grains or aggregates 
{Blum, 2004). Moreover, impact charging and successive 
charge separation can cause an electric discharge in the neb- 
ula gas. For nebula lighting {Desch and Cuzzi, 2000) a few 
hundred to thousand elementary charges per dust grain are 
required. This corresponds to impact velocities in the range 
20 . . . 100 m/s {Poppe and Schrapler, 2005) which seems 
rather high for mm particles. 

Electrostatic attraction by dipole-dipole forces has been 
seen to be important for grains of several hundred micron 
radius (chondrule size) forming clumps that are centimeters 
to tens of centimeters across {Marshall et ai, 2005; Ivlev et 
ai, 2002). Spot charges distributed over the grain surfaces 
lead to a net dipole of the grains, with growth dynamics 
very similar to that of magnetic grains. Experiments in mi- 
crogravity have shown spontaneous aggregation of particles 
in the several hundred micron size regime {Marshall and 
Cuzzi, 2001; Marshall et ai, 2005; Love and Pettit, 2004). 
The aggregates show greatly enhanced stability, consistent 
with cohesive forces increased by factors of 10^ compared 
to the normal van der Waals interaction. Based on the ex- 
periments, for weakly charged dust grains, the electrostatic 
interaction energy at contact for the charge-dipole interac- 
tion is in most cases larger than that for the charge-charge 



interaction. For heavily-charged particles, the mean mass 
of the system does not grow faster than linearly with time, 
i.e. even slower than in the non-charged case for Brownian 
motion {Ivlev etai, 2002; Konopka et ai, 2005). 

2.3. Growth and compaction of large dust aggregates 

2.3.1 Physical properties of macroscopic dust aggregates. 
Macroscopic dust aggregates can be created in the labora- 
tory by a process termed random ballistic deposition (RBD, 
Blum and Schrapler, 2004). In its idealized form, RBD uses 
individual, spherical and monodisperse grains which are de- 
posited randomly but unidirectionally on a semi-infinite tar- 
get aggregate. The volume filhng factor ^ = 0.11 of these 
aggregates, defined as the fraction of the volume filled by 
dust grains, is identical to ballistic particle-cluster aggrega- 
tion which occurs when a bimodal size distribution of par- 
ticles (aggregates of one size and individual dust grains) is 
present and when the aggregation rates between the large 
aggregates and the small particles exceed those between 
all other combinations of particle sizes. When using ide- 
alized experimental parameters, i.e. monodisperse spheri- 
cal Si02 grains with 0.75 /xm radius, Blum and Schrapler 
(2004) measured a mean volume filling factor for their 
macroscopic (cm-sized) RBD dust aggregates of (f) = 0.15, 
in full agreement with numerical predictions {Watson et 
ai 1997). Relaxing the idealized grain morphology re- 
sulted in a decrease of the volume filling factor to values of 
= 0.10 for quasi-monodisperse, irregular diamond grains 
and (/) = 0.07 for polydisperse, irregular Si02 grains {Blum, 
2004). 

Static uniaxial compression experiments with the macro- 
scopic RBD dust aggregates consisting of monodisperse 
spherical grains {Blum and Schrapler, 2004) showed that 
the volume filling factor remains constant as long as the 
stress on the sample is below ~ 500 N m~^. For higher 
stresses, the volume filling factor monotonically increases 
from (/) = 0.15 to = 0.34. Above ~ 10^ N m'^, the 
volume filUng factor remains constant at <^ = 0.33. Thus, 
the compressive strength of the uncompressed sample is 
E « 500 N m~^. These values differ from those derived 
with the models of Greenberg et ai (1995) and Sirono and 
Greenherg (2000) by a factor of a few. The compressive 
strengths of the macroscopic dust aggregates consisting of 
irregular and polydisperse grains was sUghtly lower at S 
200 N m~^. The maximum compression of these bodies 
was reached for stresses above ~ 5-10^ N and resulted 
in volume filUng factors as low as ^ = 0.20 {Blum, 2004). 
As a maximum compressive stress of ^ 10^ . . . 10^ N 
corresponds to impact velocities of ~ 15 ... 50 m/s which 
are typical for meter-sized protoplanetary dust aggregates, 
we expect a maximum volume filling factor for these bod- 
ies in the solar nebula of </> = 0.20 . . . 0.34. Blum and 
Schrapler (2004) also measured the tensile strength of their 
aggregates and found for slightly compressed samples (0 = 
0.23) T = 1,000 N m-2. Depending on the grain shape 
and the size distribution, the tensile strength decreased to 
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values of T ~ 200 N m ^ for the uncompressed case 
(Blum, 2004). 

Sirono (2004) used the above continuum properties of 
macroscopic dust aggregates, i.e. compressive strength and 
tensile strength, to model the collisions between protoplan- 
etary dust aggregates. For sticking to occur in an aggregate- 
aggregate collision, Sirono (2004) found that the impact ve- 
locity must follow the relation 



V < 0.044 



dp{(t>) 



(9) 



where p{4>) = po-4>is the mass density of the aggregate and 
po denotes the mass density of the grain material. More- 
over, the conditions 11(0) < and E((/)) < T(0) must 
be fulfilled. For the shear strength, Sirono (2004) appUes 
Y{(t>) = i,/2E((/))T((^)/3. A low compressive strength of 
the colhding aggregates favors compaction and, thus, dam- 
age restoration which can otherwise lead to a break-up of 
the aggregates. In addition, a large tensile strength also pre- 
vents the aggregates from being disrupted in the colhsion. 

Blum and Schrdpler (2004) found an approximate rela- 
tion between compressive strength and volume filUng factor 



0.8 



(10) 



which is valid in the range t^o = 0.15 < (j) < 0.21. Such a 
scaling law was also found for other types of macroscopic 
aggregates, e.g. for jammed toner particles in fluidized bed 
experiments (Valverde et al., 2004). For the aggregates con- 
sisting of monodisperse Si02 spheres, the scaling factor Eg 
can be determined to be Eg = 2.9 • N m~^. If we 
apply Eq. 10 to Eq. 9, we get, with p{(f)) = p^ ■ (p and 
Po = 2 • 10^ kg m~^, for the impact velocity of low-density 
dust aggregates 



V < 0.04. 



0.8E, 



Po{(t>-4>o) 



0.2 



0.14 {(f) - <^o)" m/s . 



(11) 

Although the function in Eq. 1 1 goes to infinity for — > 
^0. for all practical purposes the characteristic velocity is 
strongly restricted. For volume filling factors </> > 0.16 we 
get V < 0.22 m/s. Thus, following the SPH simulations by 
Sirono (2004), we expect aggregate sticking in colUsions 
for impact velocities t; < 0.2 m/s. 

2.3.2 Low-velocity collisions between macroscopic dust ag- 
gregates. Let us now consider recent results in the field of 
high-porosity aggregate colhsions. Langkowski and Blum 
(unpublished data) performed microgravity collision exper- 
iments between 0.1-1 mm-sized (projectile) RBD aggre- 
gates and 2.5 cm-sized (target) RBD aggregates. Both ag- 
gregates consisted of monodisperse spherical Si02 grains 
with radii of a = 0.75 fim. In addition to that, impact exper- 
iments with high-porosity aggregates consisting of irregular 
and/or polydisperse grains were performed. The parameter 
space of the impact experiments by Langkowski and Blum 
encompassed collision velocities in the range < v < 




Fig. 3. — (a) Non-sticking (v=1.8m/s) oblique impact between 
high-porosity dust aggregates (Langkowski and Blum, unpub- 
lished data). The three images show, from left to right, the ap- 
proaching projectile aggregate before, during, and after impact. 
The arrow in the left image denotes the impact direction of the 
projectile aggregate. The experiment was performed under micro- 
gravity conditions. It is clearly visible that the rebounding aggre- 
gate (right image) is more massive than before the collision, (b) 
Result of a high- velocity normal impact (v — 23.7 m/s) between 
compacted aggregates (Wurm et at, 2005). About half of the pro- 
jectile mass sticks to the target after the impact and is visible by 
its pyramidal structure on the flat target. Mind the different size 
scales in (a) and (b). 



3 m/s and projectile masses of 10~^ kg < m < 5 - 10~^ kg 
for all possible impact parameters (i.e. from normal to tan- 
gential impact). Surprisingly, through most of the param- 
eter space, the collisions did lead to sticking. The experi- 
ments with aggregates consisting of monodisperse spheri- 
cal Si02 grains show, however, a steep decrease in sticking 
probability from /3 = 1 to /3 = if the tangential compo- 
nent of the impact energy exceeds ~ 10~^ J (see the ex- 
ample of a non-sticking impact in Fig. 3a). Other materials 
also show the tendency towards lower sticking probabihties 
with increasing tangential impact energies. As these ag- 
gregates are "softer", the decline in sticking probability in 
the investigated parameter space is not complete. When the 
projectile aggregates did not stick to the target aggregate, 
considerable mass transfer from the target to the projectile 
aggregate takes place during the impact (Langkowski and 
Blum, unpubhshed data). Typically, the mass of the pro- 
jectile aggregate was doubled after a non-sticking colhsion 
(see Fig. 3a). 

The occurrence of sticking in aggregate-aggregate colli- 
sions at velocities > 1 m/s is clearly in disagreement with 
the prediction by Sirono (2004) (see Eq. 9). In addition, 
the evaluation of the experimental data shows that the con- 
dition for sticking, E((/)) < ¥{<}>), seems not to be fulfilled 
for high-porosity dust aggregates. This means that the con- 
tinuum aggregate model by Sirono (2004) is still not precise 
enough to fully describe the collision and sticking behavior 
of macroscopic dust aggregates. 

2.3.3 High-velocity collisions between macroscopic dust 
aggregates. The experiments described above indicate that 
at velocities above approximately 1 m/s, collisions turns 
from sticking to bouncing, at least for oblique impacts. At 
higher velocities one would naively expect that bouncing 
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and eventually erosion will continue to dominate, and this 
is also observed in a number of different experiments {Col- 
well, 2003; Bridges et al. 1996, Kouchi et al, 2002; Blum 
andMunch, 1993; Blum and Wurm, 2000). 

Growth models which assume sticking at velocities ^ 
1 m/s are therefore often considered to be impossible (e.g. 
Youdin, 2004). As velocities > 10 m/s clearly occur for 
particles that have exceeded m-size, this is a fundamental 
problem for the formation of planetesimals. 

However, recent experiments (Wurm et al., 2005) have 
studied impacts of mm-sized compact dust aggregates onto 
cm-sized compact aggregate targets at impact velocities be- 
tween 6 and 25 m/s. Compact aggregates can be the result 
of previous sticking or non-slicking collisions (see Sections 
2.3. 1 and 2.3.2). Both projectile and target consisted of /zm- 
sized dust particles. In agreement with the usual findings at 
lower impact velocities around a few m/s, the projectiles 
just rebound, shghtly fracture or even remove some parts 
of the target. However, as the velocity increases above a 
threshold of 13 m/s, about half of the mass of the projec- 
tile rigidly sticks to the target after the collision while es- 
sentially no mass is removed from the target (see Fig. 3b). 
Obviously, higher collision velocities can be favorable for 
growth, probably by destroying the internal structure of the 
porous material and dissipating energy in this way. 

Only about half of the impactor contributes to the growth 
of the target in the experiments. The other half is ejected in 
the form of small fragments, with the important implication 
that these collisions both lead to net growth of the target and 
return small particles to the disk. This keeps dust abundant 
in the disk over a long time. For the specific experiments by 
Wurm et al. (2005), the fragments were evenly distributed 
in size up to 0.5 mm. In a certain sense, the disk might thus 
quickly turn into a "debris disk" already at early times. We 
will get back to this point in section 4.4. 

3. PARTICLE-GAS INTERACTION 

Above we have seen that small solid particles grow 
rapidly into aggregates of quite substantial sizes, while re- 
taining their fractal nature (in the early growth stage) or a 
moderate to high porosity (for later growth stages). From 
the properties of primitive meteorites, we have a some- 
what different picture of nebula particulates - most of the 
soUds (chondrules, CAIs, metal grains, etc) were individu- 
ally compacted as the result of unknown melting processes, 
and were highly size-sorted. Even the porosity of what 
seem to be fine-grained accretion rims on chondrules is 25% 
or less (Scott etal., 1996; Cuzzi, 2004; Wasson etal, 2005). 
Because age-dating of chondrules and chondrites implies a 
delay of a Myr or more after formation of the first soUds, it 
seems possible that, in the asteroid formation region at least, 
widespread accretion to parent body sizes did not occur un- 
til after the mystery melting events began which formed the 
chondrules. 

It may be that conditions differed between the inner and 
outer solar system. Chondrule formation might not have 



occurred at all in the outer solar system where comet nuclei 
formed, so some evidence of the fractal aggregate growth 
stage may remain in the granular structure of comet nuclei. 
New results from Deep Impact imply that comet Tempel 
1 has a porosity of 60-80% (A'Hearn et al, 2005)! This 
value is in agreement with similar porosities found in sev- 
eral other comets {Davidsson, 2006). Even in the terrestrial 
planet/asteroid belt region, there is little reason to doubt that 
growth of aggregates started well before the chondrule for- 
mation era, and continued into and (probably) throughout 
it. Perhaps, after chondrules formed, previously ineffective 
growth processes might have dominated (sections 3.2 and 
3.3). 

3.1. Radial and vertical evolution of solids 

3.1.1 Evolution prior to formation of a dense midplane 
layer. The nebula gas (but not the particles) experiences ra- 
dial pressure gradients because of changing gas density and 
temperature. These pressure gradients act as small modifi- 
cations to the central gravity from the star that dominates 
orbital motion, so that the gas and particles orbit at differ- 
ent speeds and a gas drag force exists between them which 
constantly changes their orbital energy and angular momen- 
tum. Because the overall nebula pressure gradient force 
is outward, it counteracts a small amount of the inward 
gravitational force and the gas generally orbits more slowly 
than the particles, so the particles experience a headwind 
which saps their orbital energy, and the dominant particle 
drift is inward. Early work on gas-drag related drift was by 
Whipple (1972), Adachi et al. (1976), and Weidenschilling 
(1977). Analytical solutions for how particles interact with 
a non-turbulent nebula having a typically outward pressure 
gradient were developed by A^flfcflgflvra ef fl/. (1986). For 
instance, the ratio of the pressure gradient force to the dom- 
inant central gravity is 77 ~ 2 x 10~^, leading to a net ve- 
locity difference between the gas and particles orbiting at 
Keplerian velocity Vj^ of tjVy^ (see, e.g. Nakagawa et al. 
1986). However, if local radial maxima in gas pressure ex- 
ist, particles will drift towards their centers from both sides, 
possibly leading to radial bands of enhancement of sohds 
(see section 3.4.1). 

Small particles generally drift slowly inwards, at perhaps 
a few cm/s; even this slow inexorable drift has generated 
some concern over the years, as to how CAIs (early, high- 
temperature condensates) can survive over the apparent 1-3 
Myr period between their creation and the time they were 
incorporated into chondrite meteorite parent bodies. This 
concern, however, neglected the role of turbulent diffusion 
(see Section 3.1.2). Particles of meter size drift inwards 
very rapidly - 1 AU/century. It has often been assumed that 
these particles were "lost into the sun", but more realisti- 
cally, their inward drift first brings them into regions warm 
enough to evaporate their primary constituents, which then 
become entrained in the more slowly evolving gas and in- 
crease in relative abundance as inward migration of solids 
supplies material faster than it can be removed. Early mod- 
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els describing significant global redistribution of solids rel- 
ative to the nebula gas by radial drift were presented by 
Morfill and Volk (1984) and Stepinski and Valageas (1996, 
1997); these models either ignored midplane settUng or 
made simplifying approximations regarding it, and did not 
emphasize the potential for enhancing material in the vapor 
phase. Indeed, however, because of the large mass fluxes 
involved, this "evaporation front" effect can alter the neb- 
ula composition and chemistry significantly {Cuzzi et ai, 
2003; Cuzzi and Zahnle, 2004; Yurimoto and Kuramoto, 
2004; Krot et ai, 2005; Ciesla and Cuzzi, 2005); see also 
Cyr et al. (1999) for a discussion; however, the results of 
this paper are inconsistent with similar work by Supulver 
and Lin (2000) and Ciesla and Cuzzi (2005). This stage can 
occur very early in nebula history, long before formation of 
objects large enough to be meteorite parent bodies. 

3.1.2 The role of turbulence. The presence or absence 
of gas turbulence plays a critical role in the evolution of 
nebula sohds. There is currently no widespread agree- 
ment on just how the nebula gas may be maintained in a 
turbulent state across all regions of interest, if indeed it 
is {Stone et al, 2000, Cuzzi and Weidenschilling, 2005). 
Therefore we discuss both turbulent and non-turbulent sit- 
uations. For simplicity we will treat turbulent diffusivity T> 
as equal to turbulent viscosity i^t = acH, where c and H 
are the nebula sound speed and vertical scale height, and 
a <C 1 is a non-dimensional scaling parameter Evolution- 
ary timescales of observed protoplanetary nebulae suggest 
that 10~^ < a < 10~^ in some global sense. The largest 
eddies in turbulence have scale sizes H^/a and velocities 
Vturh = c-^a (Shakura et at., 1978; Cuzzi et at., 2001). 

Particles respond to forcing by eddies of different fre- 
quency and velocity as described by Volk et al. (1980) and 
Markiewicz et al. (1991), determining their relative veloc- 
ities with respect to the gas and to each other. The diffu- 
sive properties of MRI turbulence, at least, seem not to dif- 
fer in any significant way from the standard homogeneous, 
isotropic models in this regard {Johansen and Klahr, 2005). 
Analytical solutions for resulting particle velocities in these 
regimes were derived by Cuzzi and Hogan (2003). These 
are discussed in more detail below and by Cuzzi and Wei- 
denschilling (2005). 

Vertical turbulent diffusion at intensity a maintains par- 
ticles of stopping time in a layer of thickness h 
H^/a/Qts (Dubrulle et al, 1995; Cuzzi et al, 1996), or 
a solid density enhancement H/ h = \/Vtt^fa above the 
average value. For particles of 10 cm size and smaller and 
ot > ctmin — lG~^(a/lcm) (Cuzzi and Weidenschilling, 
2005), the resulting layer is much too large and dilute for 
collective particle effects to dominate gas motions, so radial 
drift and diffusion continue unabated. Outward radial dif- 
fusion relieves the long-standing worry about "loss into the 
sun" of small particles, such as CAIs, which are too small 
to sediment into any sort of midplane layer unless turbu- 
lence is vanishingly small (a <C ccmin), and allows some 
fraction of them to survive over 1 to several Myrs after their 



formation as indicated by meteoritic observations (Cuzzi et 
al, 2003). A similar effect might help explain the presence 
of crystalhne silicates in comets (Bockelee-Morvan, 2002; 
Gail, 2004). 

3.1.3 Dense midplane layers. When particles are able to 
settle to the midplane, the particle density gets large enough 
to dominate the motions of the local gas. This is the regime 
of collective effects; that is, the behavior of a particle de- 
pends indirectly on how all other local particles combined 
affect the gas in which they move. In regions where collec- 
tive effects are important, the mass-dominant particles can 
drive the entrained gas to orbit at nearly Keplerian velocities 
(if they are sufficiently well coupled to the gas), and thus 
the headwind the gas can exert upon the particles dimin- 
ishes from rjV^ (section 3.1.1). This causes the headwind- 
driven radial drift and all other differential particle veloci- 
ties caused by gas drag to diminish as well. 

The particle mass loading Pp/pg cannot increase with- 
out limit as particles settle, even if the global turbulence 
vanishes, and the density of settled particle layers is some- 
what self-limiting. The relative velocity solutions of Nak- 
agawa et al (1986) apply in particle-laden regimes once 
pp/Pg is known, but do not provide for a fully self consis- 
tent determination of Pp/pg in the above sense; this was ad- 
dressed by Weidenschilling (1980) and subsequently Cuzzi 
etal. (1993), Champney et al. (1995), md Dobrovolskis et 
al. (1999). The latter numerical models are similar in spirit 
to the simple analytical solutions of Dubrulle et al (1995) 
mentioned earUer, but treat large particle, high mass load- 
ing regimes in globally nonturbulent nebulae which the an- 
alytic solutions cannot address. Basically, as the midplane 
particle density increases, local, entrained gas is accelerated 
to near-Keplerian velocities by drag forces from the parti- 
cles. Well above the dense midplane, the gas still orbits at 
its pressure-supported, sub-Keplerian rate. Thus there is a 
vertical shear gradient in the orbital velocity of the gas, and 
the velocity shear creates turbulence which stirs the parti- 
cles. This is sometimes called "self-generated turbulence". 
Ultimately a steady-state condition arises where the particle 
layer thickness reaches an equilibrium between downward 
settUng and upward diffusion. This effect acts to block a 
number of gravitational instabiUty mechanisms in the mid- 
plane (section 3.3.2). 

3.2. Relative velocities and growth in turbulent and 
nonturbulent nebulae 

In both turbulent and nonturbulent regimes, particle rela- 
tive velocities drive growth to larger sizes. Below we show 
that relative velocities in both turbulent and nonturbulent 
regimes are probably small enough for accretion and growth 
to be commonplace and rapid, at least until particles reach 
meter size or so. We only present results here for particles 
up to a meter or so in size, because the expression for gas 
drag takes on a different form at larger sizes. As particles 
grow, their mass per unit area increases so they are less eas- 
ily influenced by the gas, and "decouple" from it. Their 
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overall drift velocities and relative velocities all diminish 
roughly in a linear fashion with particle radius larger than a 
meter or so (Cuzzi and Weidenschilling, 2005). 

We use particle velocities relative to the gas as derived 
by Nakagawa et al. (1986) for a range of local particle 
mass density relative to the gas density (their equations 
2.11,2.12, and 2.2 1 ) to derive particle velocities relative to 
each other in the same envirormient; all relative velocities 
scale with rjVj^. 

For simplicity we will assume particles which differ by a 
factor of three in radius; Weidenschilling (1997) finds mass 
accretion to be dominated by size spreads on this order; the 
results are insensitive to this factor. Relative velocities for 
particles of radii a and a/3, in the absence of turbulence and 
due only to differential, pressure-gradient-driven gas drag, 
are plotted in the top two panels of Fig. 4. In the top panel 
we show cases where collective effects are negligible (par- 
ticle density pp <C gas density pg). Differential vertical set- 
thng (shown at different heights z above the midplane, as 
normalized by the gas vertical scale height H) dominates 
relative velocities and particle growth high above the mid- 
plane {z/H > 0.1), and radial relative velocities dominate 
at lower elevations. Except for the largest particles, relative 
velocities for particles with this size difference are much 
less than rjVj^; particles closer in mass would have even 
smaller relative velocities. 

Moreover, in a dense midplane layer, when collective ef- 
fects dominate (section 3.1.3), all these relative velocities 
are reduced considerably from the values shown. In the 
second panel we show radial and azimuthal relative veloci- 
ties for several values of Pp/pg- When the particle density 
exceeds the gas density, collective effects reduce the head- 
wind, and all relative velocities diminish. 

Relative velocities in turbulence of several different in- 
tensities, as constrained by the nebula a (again for particles 
of radii a and a/3), are shown in the two bottom panels. 
In the second panel from the bottom, relative velocities are 
calculated as the difference of their velocities relative to the 
turbulent gas, neglecting systematic drifts and using analyti- 
cal solutions derived by Cuzzi andHogan (2003; their equa- 
tion 20) to the formalism of Volk et al. (1980). Here, the rel- 
ative velocities are forced by turbulent eddies with a range 
of size scales, having eddy turnover times ranging from the 
orbit period (for the large eddies) to much smaller values 
(for the smaller eddies), and scale with titurb = ca^/^. 

In the bottom panel we sum the various relative veloci- 
ties in quadrature to get an idea of total relative velocities 
in a turbulent nebula in which particles are also evolving by 
systematic gas-pressure-gradient driven drift. This primar- 
ily increases the relative velocities of the larger particles in 
the lower a cases. 

Overall, keeping in mind the critical velocities for stick- 
ing discussed in section 2 (^ m/s), and that particle sur- 
faces are surely crushy and dissipative, one sees that for 
particles up to a meter or so, growth by sticking is plausible 
even in turbulent nebulae for a wide range of a. Crushy ag- 
gregates will grow by accumulating smaller crushy aggre- 
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Fig. 4. — Relative velocities between particles of radii a and a/3, 
in nebulae which are non-turbulent (top two panels) or turbulent 
(bottom two panels), for a minimum mass solar nebula at 2.5AU. 
In the top panel, the particle density pp is assumed to be much 
smaller than the gas density pg: Pp/pg -C 1. Shown are the radial 
(solid line) and azimuthal (dotted line) components of the relative 
velocities. The dashed curves show the vertical relative velocities, 
which depend on height above the midplane and are shown for dif- 
ferent values of z/H. For nonturbulent cases, particles settle into 
dense midplane layers (section 3.1.3), so a more realistic situa- 
tion would be Pp/pg -c^ 1 or even ^ 1 {Cuzzi et al, 1993); thus 
in the second panel we show relative radial and angular velocities 
for three different values of Pp/pg = 0, 1, and 10. For these high 
mass loadings, z/H must be small, so the vertical velocities are 
smaller than the radial velocities. In the third panel we show rel- 
ative velocities for the same particle size difference due only to 
turbulence, for several values of a. In the bottom panel, we show 
the quadrature sum of turbulent and non-turbulent velocities, as- 
suming z/H = 0.01. 
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gates as described in earlier sections (eg. Weidenschilling, 
1997, for the laminar case). After this burst of initial growth 
to roughly meter size, however, the evolution of solids is 
very sensitive to the presence or absence of global nebula 
turbulence, as described in sections 3.3-3.4 below. Meter- 
size particles inevitably couple to the largest eddies, with 
wturb ^ several meters per second, and would destroy each 
other if they were to collide. We refer to this as the fragmen- 
tation limit. However, if particles can somehow grow their 
way past 10 meters in size, their survival becomes more as- 
sured because all relative velocities, such as shown in Fig. 4, 
decrease linearly with aps for values larger than shown in 
the plot due to the linear decrease of the area/mass ratio. 

3.2.1 Another role of gas in growth beyond the fragmen- 
tation limit. The role of gas in protoplanetary disks is 
not restricted to generate relative velocities between two 
bodies which then collide. The gas also plays an impor- 
tant role during individual collisions. A large body which 
moves through the disk faces a headwind and colhsions 
with smaller aggregates take place at its front (headwind) 
side. Fragments are thus ejected against the wind and can 
be driven back to the surface by the gas flow. 

For small bodies the gas flow can be regarded as free 
molecular flow. Thus streamlines end on the target surface 
and the gas drag is always towards the surface. Whether a 
fragment returns to the surface depends on its gas-particle 
coupling time (i.e. size and density) and on the ejection 
speed and angle. Whether reaccretion of enough fragments 
for net growth occurs, eventually depends on the distribu- 
tion of ejecta parameters, gas density, and target size. It was 
shown by Wurm et al. (2001) that growth of a larger body 
due to impact of dust aggregates entrained in a head wind 
is possible for collision velocities above 12m/s. At 1 AU a 
30-cm body in a disk model according to Weidenschilling 
and Cuzzi (1993) can grow in a coUision with small dust 
aggregates even if the initial collision is rather destructive. 

Sekiya and Takeda (2003) and KUnzli and Benz (2003) 
showed that the mechanism of aerodynamic reaccretion 
might be restricted to a maximum size due to a change in the 
flow regime from molecular to hydrodynamic. Fragments 
are then transported around the target rather then back to it. 
Wurm et al. (2004b) argue that very porous targets would 
allow some flow going through the body, which would still 
allow aerodynamic reaccretion, but this strongly depends 
on the morphology of the body (Sekiya and Takeda, 2005). 
As the gas density decreases outwards in protoplanetary 
disks, the maximum size for aerodynamic reaccretion in- 
creases. However, the minimum size also increases and the 
mechanism is only important for objects which have already 
grown beyond the fragmentation limit in some other way - 
e.g. by immediate sticking of parts of larger particles as 
discussed above (Wurm et al, 2005). 

3.3. Planetesimal formation in a midplane layer 

3.3.1 Incremental growth. Based on relative velocity argu- 
ments such as given above, Weidenschilling (1988, 1997) 



and Dullemond and Dominik (2004, 2005) find that growth 
to meter size is rapid (100-1000 yr at 1 AU; 6-7x10^ yrs 
at 30 AU) whether the nebula is turbulent or not. Such 
large particles settle towards the midplane within an orbit 
period or so. However, in turbulence, even meter-sized par- 
ticles are dispersed sufficiently that the midplane density 
remains low, and growth remains slow. A combination of 
rapid radial drift, generally erosive, high-velocity impacts 
with smaUer particles, and occasional destructive collisions 
with other meter-sized particles frustrates growth beyond 
meter-size or so under these conditions. 

In nonturbulent nebulae, even smaller particles can settle 
into fairly thin midplane layers and the total particle densi- 
ties can easily become large enough for collective effects 
to drive the entrained midplane gas to Keplerian, diminish- 
ing both headwind-induced radial drift and relative veloc- 
ities. In this situation, meter-sized particles quickly grow 
their way out of their troublesome tendency to drift radially 
(Cuzzi et al, 1993); planetesimal-sized objects form in only 
10^ - 10" years at 1 AU (Weidenschilling, 2000), and a few 
times 10^ years at 30 AU (Weidenschilling, 1997). How- 
ever, such robust growth may, in fact, be too rapid to match 
observations of several kinds (see section 4.4 and chapters 
by Dullemond et al. and Natta et al.) 

3.3.2 Particle layer instabilities. While to some workers the 
simplicity of "incremental growth" by sticking in the dense 
midplane layer of a nonturbulent nebula is appealing, past 
uncertainty in sticking properties has led others to pursue 
instability mechanisms for particle growth which are insen- 
sitive to these uncertainties. Nearly all instability mecha- 
nisms discussed to date (Safronov, 1969, 1991; Goldreich 
and Ward, 1973; Ward, 1976, 2000; SeUya, 1983, 1998; 
Goodman and Pindor, 2000; Youdin and Shu, 2002) occur 
only in nebulae where turbulence is essentially absent, and 
particle relative velocities are already very low. Just how 
low the global turbulence must be depends on the particle 
size involved, and the nebula a (sections 3.1.1 and 3.1.3). 

Classical treatments (the best known is Goldreich and 
Ward, 1973) assume that gas pressure plays no role in grav- 
itational instability, being replaced by an effective pres- 
sure due to particle random velocities (below we note this 
is not the case). Particle random velocities act to puff 
up a layer and reduce its density below the critical value, 
which is always on the order of the so-called Roche density 
p* ^ Mq /R^ where R is the distance to the central star; 
different workers give constraints which differ by factors of 
order unity (c/ Goldreich andWard, 1913, Weidenschilling, 
1980; Safronov, 1991; Cuzzi et al, 1993). These criteria 
can be traced back through Goldreich and Ward (1973) to 
Goldreich and Lynden-Bell (1965), Toomre (1964), Chan- 
drasekhar (1961) and Jeans (1928), and in parallel through 
Safronov (1960), Bel and Schatzman (1958), and Gurevitch 
and Lebedinsky (1950). Substituting typical values one de- 
rives a formal, nominal requirement that the local particle 
mass density must exceed about 10" g cm"-^ at 2 AU from 
a solar mass star even for marginal gravitational instability 
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- temporary gravitational clumping of small amplitude - to 
occur. This is about 10^ times larger than the gas density 
of typical minimum mass nebulae, requiring enhancement 
of the soUds by a factor of about 10^ for a typical average 
soUds-to-gas ratio. From section 3.1.2 we thus require the 
particle layer to have a thickness h < 10~^H, which in turn 
places constraints on the particle random velocities hfl and 
on the global value of a. 

Even assuming global turbulence to vanish, Weiden- 
schilling (1980, 1984) noted that turbulence stirred by the 
very dense particle layer itself will puff it up to thicknesses 
h that precluded even this marginal gravitational instabil- 
ity. This is because turbulent eddies induced by the vertical 
velocity profile of the gas (section 3.1.3) excite random ve- 
locities in the particles, diffusing the layer and preventing it 
from settling into a sufficiently dense state. Detailed two- 
phase fluid models by Cuzzi et al. (1993), Champney et 
al. (1995), and Dobwvolskis et al. (1999) confirmed this 
behavior. 

It is sometimes assumed that ongoing, but slow, particle 
growth to larger particles, with lower relative velocities and 
thus thinner layers (section 3.2), can lead to pp ~ p* and 
gravitational instability can then occur. However, merely 
achieving the formal requirement for marginal gravitational 
instability does not inevitably lead to planetesimals. For 
particles which are large enough to settle into suitably dense 
layers for marginal instability under self-generated turbu- 
lence {Weidenschilling, 1980; Cuzzi et al., 1993) random 
velocities are not damped on a collapse timescale, so incipi- 
ent instabilities merely "bounce" and tidally diverge. This is 
like the behavior seen in Saturn's A ring, much of which is 
gravitationally unstable by these same criteria (Salo, 1992; 
Karjalainen and Salo, 2004). Direct collapse to planetesi- 
mals is much harder to achieve, requiring much lower rel- 
ative velocities, and is unlikely to have occurred this way 
{Cuzzi et al. 1994, Weidenschilling, 1995; Cuzzi and Wei- 
denschilling, 2005). Recent results by Tanga et al. (2004) 
assume an artificial damping by gas drag and find gravi- 
tationally bound clumps form which, while not collapsing 
directly to planetesimals, retain their identity for extended 
periods, perhaps allowing for slow shrinkage; this is worth 
further numerical modeling with more reahstic damping 
physics, but still presumes a globally laminar nebula. 

For very small particles (a < 1mm; the highly relevant 
chondrule size), a different type of instability comes into 
play because the particles are firmly trapped to the gas by 
their short stopping times, and the combined system forms a 
single "one-phase" fluid which is stabilized against produc- 
ing turbulence by its vertical density gradient (Sekiya, 1998; 
Youdin and Shu, 2002; Youdin and Chiang, 2004; Garaud 
and Lin, 2004). Even for midplane layers of such small 
particles to approach a suitable density for this to occur re- 
quires nebula turbulence to drop to what may be implausi- 
bly low values (a < 10~® to 10~^°). Moreover, such one- 
phase layers, with particle stopping times tg much less than 
the dynamical collapse time (Cpp)"^/^, cannot become 
"unstable" and collapse on the dynamical timescale as nor- 



mally envisioned, because of gas pressure support, which 
is usually ignored {Sekiya, 1983; Safronov, 1991). Sekiya 
(1983) finds that particle densities must exceed lO^p* for 
such particles to undergo instabiUty and actually collapse. 
While especially difficult on one-phase instabilities by def- 
inition, this obstacle should be considered for any particle 
with stopping time much shorter than the dynamical col- 
lapse time - that is, pretty much anything smaller than a 
meter for pp p*. 

A slower "sedimentation" from axisymmetric rings (or 
even localized blobs of high density, which might form 
through fragmentation of such dense, differentially rotat- 
ing rings), has also been proposed to occur under condi- 
tions normally ascribed to marginal gravitational instability 
{Sekiya, 1983; Safronov, 1991; Ward, 2000), but this ef- 
fect has only been modeled under nonturbulent conditions 
where, as mentioned above, growth can be quite fast by 
sticking alone. In a turbulent nebula, diffusion (or other 
complications discussed below, such as large vortices, spi- 
ral density waves, etc) might preclude formation of all but 
the broadest- scale "rings" of this sort, which have radial 
scales comparable to H and grow only on extremely long 
timescales. 

3.4. Planetesimal formation in turbulence 

A case can be made that astronomical, asteroidal, and 
meteoritic observations require planetesimal growth to stall 
at sizes much smaller than several km, for something Uke 
a million years {DuUemond and Dominik, 2005; Cuzzi and 
Weidenschilling, 2005; Cuzzi et al, 2005). This is perhaps 
most easily explained by the presence of ubiquitous weak 
turbulence {a > 10"'*). Once having grown to meter-size, 
particles couple to the largest, most energetic turbulent ed- 
dies, leading to mutual colUsions at relative velocities on 
the order of uturb ~ \/ctc ^ 30 m/s, which are proba- 
bly disruptive, staUing incremental growth by sticking at 
around a meter in size. Astrophysical observations support- 
ing this inference are discussed in the next section. In prin- 
ciple, planetesimal formation could merely await cessation 
of nebula turbulence and then happen all at once; pros and 
cons of this simple concept are discussed by Cuzzi and Wei- 
denschilling (2005). The main difficulty with this concept is 
the very robust nature of growth in dense midplane layers of 
nonturbulent nebulae, compared to the very extended dura- 
tion of 10® years which apparently characterized meteorite 
parent body formation (chapter by Wadhwa et al). Further- 
more, if turbulence merely ceased at the appropriate time 
for parent body formation to begin, particles of all sizes 
would settle and accrete together, leaving unexplained the 
very well characterized chondrite size distributions we ob- 
serve. Alternately, several suggestions have been advanced 
as to how the meter- sized barrier might be overcome even 
in ongoing turbulence, as described below. 

3.4.1 Concentration of boulders in large nebula gas struc- 
tures. The speedy inward radial drift of meter-sized par- 
ticles in nebulae where settUng is precluded by turbulence 
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might be slowed if they can be, even temporarily, trapped by 
one of several possible fluid dynamical effects. It has been 
proposed that such trapping concentrates them and leads to 
planetesimal growth as well. 

Large nebula gas dynamical structures such as systemat- 
ically rotating vortices (not true turbulent eddies) have the 
property of concentrating large boulders near their centers 
(Barge and Sommeria, 1995; Tanga et al., 1996; Bracco et 
al, 1998; GodonandLivio, 2000; Klahr and Bodenheimer, 
2006). In some of these models the vortices are simply 
prescribed and/or there is no feedback from the particles. 
Moreover, there are strong vertical velocities present in re- 
aUstic vortices, and the vortical flows which concentrate m- 
size particles are not found near the midplane, where the 
m-sized particles reside (Barranco and Marcus, 2005). Fi- 
nally, there may be a tendency of particle concentrations 
formed in modeled vortices to drift out of them and/or de- 
stroy the vortex (Johansen et al., 2004). 

Another possibility of interest is the buildup of solids 
near the peaks of nearly axisymmetric, localized radial pres- 
sure maxima, which might for instance be associated with 
spiral density waves (Haghighipour and Boss, 2003a,b; 
Rice et al, 2004). Johansen et al. (2006) noted boul- 
der concentration in radial high pressure zones of their full 
simulation, but (in contrast to above suggestions about vor- 
tices), saw no concentration of meter-sized particles in the 
closest thing they could resolve in the nature of actual tur- 
bulent eddies. Perhaps this merely highlights the key differ- 
ence between systematically rotating (and often artificially 
imposed) vortical fluid structures, and realistic eddies in re- 
alistic turbulence. 

Overall, models of boulder concentration in large-scale 
fluid structures will need to assess the tendency for rapidly 
colUding meter-sized particles in such regions to destroy 
each other, in the real turbulence which will surely accom- 
pany such structures. For instance, breaking spiral density 
waves are themselves potent drivers for strong turbulence 
(Bo/ey era/., 2005). 

3.4.2 Concentration of chondrules in 3D turbulence. An- 
other suggestion for particle growth beyond a meter in 
turbulent nebulae is motivated by observed size-sorting in 
chondrites. Cuzzi et al. (1996, 2001) have advanced the 
model of turbulent concentration of chondrule- sized (mm 
or smaller diameter) particles into dense zones, that ulti- 
mately become the planetesimals we observe. This effect, 
which occurs in genuine, 3D turbulence (both in numeri- 
cal models and laboratory experiments), naturally satisfies 
meteoritics observations in several ways under quite plausi- 
ble nebula conditions. It offers the potential to leapfrog the 
problematic meter-size range entirely and would be appU- 
cable (to differing particle types) throughout the solar sys- 
tem (see Cuzzi and Weidenschilling, 2005 and Cuzzi et al., 
2005 for reviews). This scenario faces the obstacle that the 
dense, particle -rich zones which certainly do form are far 
from solid density, and might be disrupted by gas pressure 
or turbulence before they can form soUd planetesimals. As 



with dense midplane layers of small particles, gas pressure 
is a formidable barrier to gravitational instability on a dy- 
namical timescale in dense zones of chondrule- sized par- 
ticles formed by turbulent concentration. However, as with 
other small-particle scenarios, sedimentation is a possibility 
on longer timescales than that of dynamical collapse. It is 
promising that Sekiya (1983) found that zones of these den- 
sities, while "incompressible" on the dynamical timescale, 
form stable modes. Current studies are assessing whether 
the dense zones can survive perturbations long enough to 
evolve into planetesimals. 

3.5. Summary of the situation regarding planetesimal 
formation 

As of the writing of this chapter, the path to planetesimal 
formation remains unclear. In nonturbulent nebulae, a va- 
riety of options seem to exist for growth which - while not 
on dynamical collapse timescales, is rapid on cosmogonic 
timescales (^ 10^ years). However, this set of conditions 
and growth timescales seems to be at odds with asteroidal, 
meteoritic, and astronomical observations of several kinds 
(Russell etal, 2006; Dullemond and Dominik, 2005; Cuzzi 
and Weidenschilling, 2006; Cuzzi et al., 2005; chapter by 
Wadhwa et al.). The alternate set of scenarios - growth be- 
yond a meter or so in size in turbulent nebulae - are per- 
haps more consistent with the observations but are still in- 
completely developed beyond some promising directions. 
The chaUenge is to describe quantitatively the rate at which 
planetesimals form under these inefficient conditions. 

4. GLOBAL DISK MODELS WITH SETTLING 
AND AGGREGATION 

Globally modeling a protoplanetary disk including dust 
setthng, aggregation, radial drift and mixing, along with ra- 
diative transfer solutions for the disk temperature and spec- 
trum form a major numerical challenge, because of the 
many orders of magnitude that have to be covered both in 
time scales (inner disk versus outer disk, growth of small 
particles versus growth of large objects) and particle sizes. 
Further numerical difficulties result from the fact that small 
particles may contribute significantly to the growth of larger 
bodies, and careful renormalization schemes are necessary 
to treat these processes correctly and in a mass-conserving 
manner (Dullemond and Dominik, 2005). Further difficul- 
ties arise from uncertainty about the strength and spatial ex- 
tent of turbulence during the different evolutionary phases 
of a disk. A complete model covering an entire disk and the 
entire growth process along with all relevant disk physics 
is currently still out of reach. Work so far has therefore ei- 
ther focused on specific locations in the disk, or has used 
parametrized descriptions of turbulence with limited sets of 
physical growth processes. However, these "single slice" 
models have the problem that radial drift can become so 
large for m-sized objects, that these leave the slice on a 
time scale of a few orbital times (Weidenschilling, 1977; 
section 3.1.1). Nevertheless, important results have come 
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forth from these efforts, that test underlying assumptions of 
the models. 

For the spectral and imaging appearance of disks, there 
are two main processes that should produce easily observ- 
able results: particle settling and particle growth. Particle 
settling is due to the vertical component of gravity acting 
in the disk on the pressure-less dust component (section 
3.1.2). Neglecting growth for the moment, settling leads 
to a vertical stratification and size sorting in the disk. Small 
particles settle slowly and should be present in the disk at- 
mosphere for a long time, while large particles settle faster 
and to smaller scale heights. While in a laminar nebula this 
is a purely time dependent phenomenon, this result is per- 
manent in a turbulent nebula as each particle size is spread 
over its equilibrium scale height (Dubrulle et al., 1995). 
From a pure settling model, one would therefore expect that 
small dust grains will increasingly dominate dust emission 
features (cause strong feature-to-continuum ratios) as large 
grains disappear from the surface layers. 

Grain growth may have the opposite effect. While verti- 
cal mixing and settUng still should lead to a size stratifica- 
tion, particle growth can become so efficient that all small 
particles are removed from the gas. In this case, dust emis- 
sion features should be characteristic for larger particles 
(i.e. no or weak features, see chapter by Natta et al.). At the 
same time, the overall opacity decreases dramatically. This 
effect can become significant, as has been realized already 
early on (Weidenschilling, 1980, 1984; Mizuno, 1989). In 
order to keep the small particle abundance at realistic levels 
and the dust opacity high, Mizuno et al. (1988) considered a 
steady inflow of small particles into a disk. However, disks 
with signs of small particles are still observed around stars 
that seem to have completely removed their parental clouds, 
so this is not a general solution for this problem. In the fol- 
lowing we discuss the different disk models documented in 
the hterature. We begin with a discussion of earlier models 
focusing on specific regions of the solar system. 

4.1. Models limited to specific regions in the solar sys- 
tem 

Models considering dust settling and growth in a sin- 
gle vertical sUce have a long tradition, and have been re- 
viewed in previous Protostars and Planets HI {Cuzzi and 
Weidenschilling, 1993). We therefore refrain from an in- 
depth coverage and only recaU a few of the main results. 
The global models discussed later are basically similar cal- 
culations, with higher resolution, and for a large set of radii. 

Weidenschilling has studied the aggregation in lami- 
nar {Weidenschilling, 1980, 2000) and turbulent {Weiden- 
schilling, 1984, 1988) nebulae, focusing on the region 
of terrestrial planet formation, in particular around lAU. 
These papers contain the basic descriptions of dust settling 
and growth under laminar and turbulent conditions. They 
show the occurrence of a rain-out after particles have grown 
to sizes where the settUng motion starts to exceed the ther- 
mal motions. Nakagawa et al. (1981, 1986) study settling 



and growth in vertical slices, also concentrating on the ter- 
restrial planet formation regions. They find that within 3000 
years, the midplane is populated by cm-sized grains. Wei- 
denschilling (1997) studied the formation of comets in the 
outer solar system with a detailed model of a non-turbulent 
nebula, solving the coagulation equation around 30AU. In 
these calculations, growth initially proceeds by Brownian 
motion, without significant settling, for the first lOOOOyrs. 
Then, particles become large enough and start to settle, 
so that the concentration of soUds increases quickly after 
5 X lO^yr The particle layer reaches the critical density 
where the layer gravitational instabihty is often assumed to 
occur, but first the high velocity dispersion prevents the col- 
lapse. Later, a transient density enhancement still occurs, 
but due to the small coUisional cross section of the typically 
Im-sized bodies, growth must still happen in individual 2- 
body coUisions. 

4.2. Dust aggregation during early disk evolution 

Schmitt et al. (1997) implemented dust coagulation in an 
a disk model. They considered the growth of PCA in a one- 
dimensional disk model, i.e. without resolving the vertical 
structure of the disk. The evolution of the dust size distribu- 
tion is foUowed for 100 years only. In this time, at a radius 
of 30AU from the star, first the smallest particles disappear 
within 10 years, due to Brownian motion aggregation. This 
is followed by a self-similar growth phase during which the 
volume of the particles increases by 6 orders of magnitude. 
Aggregation is faster in the inner disk, and the decrease in 
opacity foUowed by rapid cooUng leads to a thermal gap 
in the disk around 3AU. Using the CCA particles, aggre- 
gation stops in this model after the small grains have been 
removed. For such particles, longer timescales are required 
to continue the growth. 

Global models of dust aggregation during the prestel- 
lar coUapse stage and into the early disk formation stage 
are numerically feasible because the growth of particles 
is limited. Suttner et al. (1999) and Suttner and Yorke 
(2001) study the evolution of dust particles in protostel- 
lar envelopes, during collapse, and the first 10* years of 
dynamical disk evolution, respectively. These very ambi- 
tious models include a radiation hydrodynamic code that 
can treat dust aggregation and shattering using an implicit 
numerical scheme. They find that during a collapse phase of 
10^ years, dust particles grow due to Brownian motion and 
differential radiative forces, and can be shattered by high 
velocity collisions cause by radiative forces. During early 
disk evolution, they find that at 30AU from the star within 
the first pressure scale height from the midplane, small par- 
ticles are heavily depleted because the high densities lead 
to frequent colhsions. The largest particles grow by a factor 
of 100 in mass. Similar results are found for PCA particles, 
while CCA particles show accelerated aggregation because 
of the enhanced cross section in massive particles. Within 
10* years, most dust moves to the size grid limit of 0.2mm. 
While aggregation is significant near the midplane (opaci- 
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ties are reduced by more than a factor 10), the overall struc- 
ture of the model is not yet affected strongly, because at the 
low densities far from the midplane aggregation is Umited 
and changes in the opacity are only due to differential ad- 
vection. 

4.3. Global settling models 

Settling of dust without growth goes much slower than 
settling that is accelerated by growth. However, even pure 
settling calculations show significant influence on the spec- 
tral energy distributions of disks. While the vertical optical 
depth is unaffected by setthng alone, the height at which 
stellar light is intercepted by the disk surface changes. 
Miyake and Nakagawa (1995) computed the effects of dust 
settling on the global SED and compared these results with 
IRAS observations. They assume that after the initial set- 
tling and growth phase, enough small particles are left in 
the disk to provide an optically thick surface and follow the 
decrease of the height of this surface, concluding that this is 
consistent with the life-times of T Tauri disks, because the 
settling time of a 0. 1 fim grain within a single pressure scale 
height is of order 10 Myr. However, the initial settUng phase 
does lead to strong effects on the SED, because settling 
times at several pressure scale heights are much shorter. 

Dullemond and Dominik (2004) show that settling from 
a fully mixed passive disk leads to a decrease of the surface 
height in 10^-10^ years, and can even lead to self-shadowed 
disks (see chapter by Dullemond et al). 

4.4. Global models of dust growth 

Mizuno (1989) computes global models including evap- 
oration, and a steady state assumption using small grains 
continuously raining down from the ISM. The vertical disk 
structure is not resolved, only a single zone in the midplane 
is considered. He finds that the Rosseland mean opacity de- 
creases, but then stays steady due to the second generation 
grains. 

Kornet et al. (2001) model the global gas and dust disk 
by assuming that at a given radius, the size distribution of 
dust particles (or planetesimals) is exactly monodisperse, 
avoiding the numerical complications of a full solution of 
the Smoluchowski equation. They find that the distribution 
of solids in the disk after 10^ years depends strongly on the 
initial mass and angular momentum of the disk. 

Ciesla and Cuzzi (2005) model the global disk using a 
four-component model: Dust grains, m-size boulders, plan- 
etesimals and the disk gas. This model tries to capture the 
main processes happening in a disk: growth of dust grains 
to m-sized bodies, the migration of m-sized bodies and the 
resulting creation of evaporation fronts, and the mixing of 
small particles and gas by turbulence. The paper focuses 
on the distribution of water in the disk, and the dust growth 
processes are handled by assuming timescales for the con- 
version from one size to the next. Such models are therefore 
mainly useful for the chemical evolution of the nebula and 
need detailed aggregation calculations as input. 
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Fig. 5. — Time evolution of the disk SED in the laminar case (top 
panel) and the turbulent case (bottom panel). From DD05 

The most complete long-term integrations of the equa- 
tions for dust settiing and growth are described in recent pa- 
pers by Tanaka et al. (2005, henceforth THI05) and Dulle- 
mond and Dominik (2005, henceforth DD05). These pa- 
pers implement dust settling and aggregation in individual 
vertical slices through a disk, and then use many slices to 
stitch together an entire disk model, with predictions for 
the resulting optical depth and SED from the developing 
disk. Both models have different Umitations. THI05 con- 
sider only laminar disk models, so that turbulent mixing 
and collisions between particles driven by turbulence are 
not considered. Their calculations are limited to compact 
solid particles. DD05's model is incomplete in that it does 
not consider the contributions of radial drift and differential 
angular velocities between different particles. But in addi- 
tion to calculations for a laminar nebula, they also introduce 
turbulent mixing and turbulent coagulation, as well as PCA 
and CCA properties for the resulting dust particles. THI05 
use a two-layer approximation for the radiative transfer so- 
lution, while DD05 run a 3D Monte-Carlo radiative transfer 
code to compute the emerging spectrum of the disk. Both 
models find that aggregation proceeds more rapidly in the 
inner regions of the disk than in the outer regions, quickly 
leading to a region of low optical depth in the inner disk. 
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Both calculations find that a bi-modal size distribution 
is formed, with large particles in the midplane, formed by 
rainout (the fast setthng of particles after their settling time 
has decreased below their growth time) and continuing to 
grow quickly, and smaller particles remaining higher up in 
the disk and then slowly trickling down. In the laminar disk, 
growth stops in the DD05 calculations at cm sizes because 
radial drift was ignored, in THI05, particles continue to 
grow beyond this regime. 

The settUng of dust causes the surface height of the disk 
to decrease, reducing the overall capacity of the disk to re- 
process stellar radiation. THI05 find that at 8 AU from the 
star the optical depth of the disk at 10/xm reaches about 
unity after a bit less than lO^yrs. In the inner disk, the sur- 
face height decreases to almost zero in less than lO^yrs. 
The SED of the model shows first a strong decrease at 
wavelength of 100/im and longer, within the first 10''yrs. 
After that the near-IR and mid-IR radiation also decreases 
sharply. THI05 consider their results to be roughly consis- 
tent wilh Ihe observations of decreasing fluxes at near-IR 
and mm wavelengths in disks. 

The calculations by DD05 show a more dramatic effect, 
as shown in Fig. 5. In the calculations for a laminar disk, 
here the surface height already significantly decreases in 
the first 10* years, then the effect on the SED is initially 
strongest in the mid-lR region. After lO^yrs, the fluxes have 
dropped globally by at least a factor of 10, except for the 
mm regime, which is affected greatly only after a few times 
10® years. 

In the calculations for a turbulent disk (DD05), the de- 
pletion of small grains in the inner disk is strongly en- 
hanced. This result is caused by several effects. First, tur- 
bulent mixing keeps the particles moving even after they 
have settled to the mid-plane, allowing them to be mixed 
up and rain down again through a cloud of particles. Fur- 
thermore, vertical mixing in the higher disk regions mixes 
low-density material down to higher densities, where aggre- 
gation can proceed much faster. The material being mixed 
back up above the disk is then largely deprived of sohds, be- 
cause the large dust particles decouple from the gas and stay 
behind, settling down to the mid plane. The changes to the 
SED caused by coagulation and settling in a turbulent disk 
are dramatic, and clearly inconsistent with the observations 
of disks around T Tauri stars that indicate lifetimes of up to 
10*^ years. DD05 conclude that ongoing particle destruction 
must play an important role, leading to a steady-state size 
distribution for small particles (section 2.3.1). 

4.5. The role of aggregate structure 

Up to now, most solutions for the aggregation equation 
in disks are still based on the assumption of compact par- 
ticles resulting from the growth process. However, at least 
for the small aggregates formed initially, this assumption 
is certainly false. First of all, if aggregates are fluffy, with 
large surface-to-mass ratios, it will be much easier to keep 
these particles in the disk surface where they can be ob- 



served as scattering and IR emitting grains. Observations 
of the 10;um silicate features show that in many disks, the 
population emitting in this wavelength range is dominated 
by particles larger than interstellar {van Boekel et al, 2005; 
Kessler-Silacci et al, 2005). When modelled with compact 
grains, the typical size of such grains is several microns, 
with corresponding settling limes less than a Myr. When 
modelled with aggregates, particles have to be much larger 
to produce similar signatures (flattened feature shapes, e.g. 
Minetai, 2005). 

When considering the growth time scales, in particular in 
regions where settling is driving the relative velocities, the 
timescales are surprisingly similar to the case of compact 
particles {Safronov, 1969, WeidenschilUng, 1980). While 
initially, fluffy particles settle and grow slowly because of 
small settling velocities, the larger coUisional cross section 
soon leads to fast collection of small particles, and fluffy 
particles reach the mid-plane as fast as compact grains, and 
with similar masses collected. 

5. SUMMARY AND FUTURE PROSPECTS 

A lot has been achieved in the last few years, and our 
understanding of dust growth has advanced significantly. 
There are a number of issues where we now have clear an- 
swers. However, a number of major controversies remain, 
and future work will be needed to address these before we 
can come to a global picture of how dust growth in pro- 
toplanetary disks proceeds and which of the possible ways 
toward planetesimals are actually used by nature. In table 1 
on the following page we summarize our main conclusions 
and questions, and note some priorities for research in the 
near future. 
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Table 1: Overview 



What We ReaUy Know 


Main Controversies/Questions 


Future Priorities 


Microphysics 


— Dust particles stick in collisions with 
less than ~ 1 m/s velocity due to van der 
Waals force or hydrogen bonding. 

— For low relative velocities (<C Im/s) a 
cloud of dust particles evolves into frac- 
tal aggregates (Dj < 2) with a quasi- 
monodisperse mass distribution. 

— Due to the increasing collision energy, 
growing fractal aggregates can no longer 
keep their structures so that non-fractal 
(but very porous) aggregates form (still at 
V <C Im/s). 

— Macroscopic aggregates have porosities 

> 65% when collisional compaction, and 
not sintering or melting occurs. 


— At what aggregate size does com- 
paction happen in a nebula environment? 

— When do collisions between macro- 
scopic aggregates result in sticking? Some 
experiments show no sticking at rather 
low impact velocities, while others show 
sticking at high impact speeds. 

— How important are special material 
properties: organics, ices, magnetic and 
electrically charged particles? 

— What are the main physical parame- 
ters (e.g. velocity, impact angle, aggre- 
gate porosity/material/shape/mass) deter- 
mining the outcome of a collision? 


— More empirical studies in collisions be- 
tween macroscopic aggregates required. 

— Macroscopic model for aggregate col- 
lisions (continuimi description) based on 
microscopic model and experimental re- 
sults. 

— Develop recipes for using the micro- 
physics in large scale aggregation calcu- 
lations. 

— Develop aggregation models that treat 
aggregate structure as a variable in a self- 
consistent way. 


Nebula processes 


— Particle velocities and relative veloci- 
ties in turbulent and nonturbulent nebulae 
are understood; values are < 1 m/s for 
ap < 1 — 3 g cm~^ depending on alpha. 

— Radial drift decouples large amounts 
of solids from the gas and migrates it ra- 
dially, changing nebula mass distribution 
and chemistry 

— Turbulent diffusion can offset inward 
drift for particles of cm size and smaller, 
relieving the "problem" about age differ- 
ences between CAIs and chondrules. 


— What happens to dust aggregates in 
highly mass-loaded regions in the solar 
nebula, e.g. midplane, eddies, stagnation 
points? 

— Is the nebula turbulent? If so, how 
does the intensity vary with location and 
time? Can purely hydrodynamical pro- 
cesses produce self-sustaining turbulence 
in the terrestrial planet formation zone? 

— Can large-scale structures (vortices, 
spiral density waves) remain stable long 
enough to concentrate boulder-size parti- 
cles? 

— Can dense turbulently concentrated 
zones of chondrule-size particles survive 

to become actual planetesimals? 


— Relative velocities in highly mass- 
loaded regions in the solar nebula, e.g. 
midplane, eddies, stagnation points. 

— Improve our understanding of turbu- 
lence production processes at very high 
nebula Reynolds numbers. 

— Model effects of MRI-active upper lay- 
ers on dense, non-ionized gas in magneti- 
cally dead zones. 

— Model the evolution of dense strength- 
less clumps of particles in turbulent gas. 

— Model collisional processes in boulder- 
rich vortices and high-pressure zones. 

— Model evolution of dense clumps in tur- 
bulent gas. 


Global modelling and comparison with observations 


— Small grains are quickly depleted by in- 
corporation into larger grains. 

— Growth timescales are short for small 
compact and fractal grains alike. 

— Vertical mixing and small grain re- 
plenishment are necessary to keep the ob- 
served disk structures (thick/flaring). 


— What is the role of fragmentation for 
the small grain component? 

— Are the "small" grains seen really large, 
fluffy aggregates? 

— Are the mm/cm sized grains seen in 
observations compact particles, or much 
larger fractal aggregates? 

— What is the global role of radial trans- 
port? 


— Study the optical properties of large ag- 
gregates, fluffy and compact. 

— Implement realistic opacities in disk 
models to produce predictions and com- 
pare with observations. 

— Construct truly global models including 
radial transport. 

— More resolved disk images at many 
wavelengths, to better constrain models. 
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